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ABSTRACT 

We model the spectral energy distribution of the type 1 Seyfert galaxy NGC 5548, fitting data from simultaneous 
optical, UV, and X-ray monitoring observations. We assume a geometry consisting of a hot central Comptonizing 
region surrounded by a thin accretion disk. The properties of the disk and the hot central region are determined by 
the feedback occurring between the hot Comptonizing region and thermal reprocessing in the disk that, along with 
viscous dissipation, provides the seed photons for the Comptonization process. The constraints imposed upon this 
model by the multiwavelength data allow us to derive limits on the central black hole mass, M <2xW M 0 , the 
accretion rate, M < 2.5 x 10^ MQyr“'/M, and the radius of the transition region between the thin outer disk and 
the geometrically thick, hot inner region, ~ 2-5 x lO'"^ cm. 

Subject headings: galaxies: active — galaxies: individual (NGC 5548) — galaxies: Seyfert — X-rays: galaxies 


1. INTRODUCTION 

In two previous papers, we developed a model of thermal 
Comptonization and disk thermal reprocessing to explain the 
data from recent multiwavelength monitoring observations of 
the type 1 Seyfert galaxies NGC 3516 and NGC 7469 (Chiang 
& Blaes 2001, hereafter Paper I; Chiang 2002, hereafter Pa¬ 
per II). The essential features of this model were proposed by 
Poutanen, Krolik, & Ryde (1997) to account for the different 
spectral states of the Galactic X-ray binary Cygnus X-1; and a 
simplified version was used by Zdziarski, Lubihski, & Smith 
(1999) to model the correlation seen in X-ray binaries and type 
1 Seyfert galaxies between the X-ray spectral index and the rel¬ 
ative strength of the Compton reflection continuum. A key el¬ 
ement of this model is the feedback occurring between a cen¬ 
trally located Comptonizing plasma and the thin accretion disk 
that encircles it. The disk reprocesses X-rays from the central 
region, thus providing a significant fraction of the seed photons 
for the Comptonization process itself. Eor a given set of opti¬ 
cal, UV, and X-ray data, this feedback mechanism determines 
the equilibrium state of the Comptonizing plasma. By comput¬ 
ing the equilibria required by the NGC 3516 and NGC 7469 
data, we found that the seemingly uncorrelated optical/UV and 
X-ray continuum variations exhibited by these sources could be 
described almost entirely in the context of this model. Eurther- 
more, in modeling these variations, we were able to constrain 
the central black hole masses and the accretion rates of these 
objects; and we were able to compute the X-ray continuum at 
energies extending beyond ~ 20 keV, a part of the spectrum 
that will soon be made more accessible by the International 
Gamma-Ray Astrophysics Laboratory (INTEGRAL). 

During an intensive three day monitoring campaign in 1998 
April with the Rossi X-ray Timing Explorer (RXTE), the Ad¬ 
vanced Satellite for Cosmology and Astrophysics (ASCA), and 
the Hubble Space Telescope, NGC 3516 exhibited an optical 
continuum that was surprisingly constant, consisting of only ^ 
3% variations, while simultaneously showing 60% changes 


in the 2-10 keV X-ray flux (Edelson et al. 2000). This appar¬ 
ent lack of response by the optical continuum to the changes in 
the X-ray flux has been taken as evidence against the thermal 
reprocessing model for the optical and UV variability. How¬ 
ever, despite the relatively small amplitude of its fluctuations, 
the optical continuum varied essentially simultaneously across 
a wide range of wavelengths, in a manner consistent with disk 
thermal reprocessing (e.g., Krolik et al. 1991; Courvoisier & 
Clavel 1991; Collin-Souffrin 1991). In Paper I, we showed 
that a constant optical continuum in the presence of large X- 
ray fluctuations could be reconciled with the disk thermal re¬ 
processing model if the radial extent, r^, of the central Comp¬ 
tonizing plasma varied inversely with the X-ray luminosity, Lj^, 
so that the flux incident upon the disk at moderately large radii 
remained constant. Slight departures from the L^ oc rela¬ 
tionship would then produce the highly correlated, low ampli¬ 
tude variations seen in the optical continuum of NGC 3516. We 
found that the range of X-ray fluxes and spectral indexes exhib¬ 
ited by this source is naturally accommodated by this behavior. 

The interpretation of the NGC 7469 data is more complex. 
Data from a 30 day monitoring campaign in 1996 June-July 
with the International Ultraviolet Explorer (lUE) and RXTE 
showed that the UV continuum at A1350A and the 2-10 keV 
X-ray light curves were not correlated in any simple way (Nan- 
dra et al. 1998). Although the amplitude and time scales of the 
most prominent variations in the two datasets were comparable, 
some features in the UV occurred nearly simultaneously with 
similar looking features in the X-rays, while other features that 
were similar in the two light curves were offset by ^ 4 days. 
A quantitative analysis revealed that the strongest correlation 
between the two light curves corresponded to the UV varia¬ 
tions leading those of the X-rays by several days, contrary to 
naive expectations of the thermal reprocessing model. Nandra 
et al. ( 2000 ) were able to explain these data, at least qualita¬ 
tively, by noting that the UV fluxes were highly correlated at 
near zero temporal lag with the X-ray spectral index and also 
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with a broad band (0.1-100 keV) extrapolation of the 2-10 keV 
X-ray power-law. When we applied the thermal Comptoniza- 
tion/disk reprocessing model to these data, we found that rather 
large changes in the bolometric X-ray luminosity were required 
in order to produce the relatively large variations seen in the 
UV (Paper II). Since the X-ray continua for this source are rel¬ 
atively soft and are constrained in the 2-10 keV band by the 
RXTEIVCA observations, a broad range of temperatures for the 
Comptonizing electrons is needed so that the thermal Comp- 
tonization process can produce a concomitantly large range of 
X-ray luminosities. 

Despite the difficulties associated with explaining these ob¬ 
servations, the application of the model to these data yielded 
some encouraging results. First, the relatively large inner trun¬ 
cation radius of the thin disk, r^in ~ 3 x 10'“* cm, that we found 
by fitting the optical and UV continuum data for NGC 7469 
is consistent with that inferred from ASCA observations of the 
narrow iron Ka emission line by Guainazzi et al. (1994) if the 
broad line region reverberation mapping estimates of the black 
hole mass are correct (Peterson & Wandel 2000). Second, the 
value we inferred for the product of the black hole mass and ac¬ 
cretion rate, MM ~ 10® Mq yr“', agrees with that found by Col¬ 
lier et al. (1998) who examined time delays between UV and 
optical continuum bands under the assumption of thermal re¬ 
processing. Third, a global energy balance condition is satisfied 
by our models for both NGC 3516 and NGC 7469 wherein the 
total bolometric luminosity is consistent with the inferred ac¬ 
cretion rate assuming nominal black hole masses of IO^Mq 
for each source. 

In the present paper, we analyze simultaneous optical, UV, 
and X-ray data for the Seyfert 1 galaxy NGC 5548 which 
were taken during two separate monitoring observations: a 
ground-based-optical/ZUE/G/nga campaign in 1989-90 (Peter¬ 
son et al. 1991; Clavel et al. 1991, 1992; Nandra et al. 1991) 
and our own 1998 ground-based-optical/Tixfreme Ultraviolet 
Explorer (EUVE)IASCA/RXTE monitoring campaign (Chiang 
et al. 2000; Dietrich et al. 2001). Using our model, these data 
allow us to place limits on the mass and accretion rate of this 
object; and from that information, we can make predictions 
for what should be seen by XMM-Newton observations of this 
source in the optical, UV, and X-rays and by possible forth¬ 
coming simultaneous observations by INTEGRAL with other 
X-ray satellites such as XMM-Newton, the Chandra X-ray Ob¬ 
servatory, and RXTE. In §2 of this paper, we describe the op¬ 
tical, UV, and X-ray data. In §3, we present our semi-analytic 
model of the spectral energy distributions (SEDs); and in §4, 
we derive constraints on the black hole mass and accretion rate. 
In §5, we show results from detailed Monte Carlo simulations 
of the thermal Comptonization and disk thermal reprocessing 
physics that verify the empirical relations we used in §3 and the 
SEDs we computed in §4. In §6, we discuss our results and 
their implications. We conclude in §7. 

2. OBSERVATIONS OF NGC 5548 

The simultaneous giound-hased/IUE/Ginga observations of 
NGC 5548 cover a period Jan 1989 through July 1990. The 
lUE observations themselves took place from 1988 Dec to 1990 
July and are described in Clavel et al. (1991) and Clavel et al. 
(1992). These references provide the continuum fluxes at the 
source rest frame wavelengths A1350A, A1840A, and A2670A. 

^ Magdziarz et al. (1998) were partly guided by non-simultaneous ROSAT data, 
residuals above a power law at the softest X-ray energies. 


Eor the fluxes at the latter two wavelengths, we make the cor¬ 
rections for Eell and Balmer continuum contamination as de¬ 
scribed by Clavel et al. (1992; see also Wamsteker et al. 1990). 
The optical fluxes at A5100A were obtained via ground-based 
observations by Peterson et al. (1991). We correct for extinc¬ 
tion for all these wave bands and for starlight contamination in 
the optical band as described in Magdziarz et al. (1998). 

The Ginga data are from the monitoring campaign reported 
by Nandra et al. (1991) and were re-extracted by Magdziarz et 
al. (1998) to obtain additional effective area at the higher ener¬ 
gies (> 15 keV). The re-extracted data cover energies from 1.7 
to 30 keV. We have fit the Ginga data from the nine epochs that 
were simultaneous with the lUE observations using a power- 
law + reflection model for the continuum components (the 
PEXRAV model [Magdziarz & Zdziarski 1995] in XSPEC), a 
Gaussian function for the fluorescent iron Ka emission line, 
and a Galactic absorbing column fixed at Nh = 1.65 x 10^® cm. 
In addition, we included an ionized absorber that is illuminated 
by a continuum with the same spectral index as the cut-off 
power-law (the ABSORI model [Done et al. 1992] in XSPEC). 
This latter component is required to model the absorption at 
energies below ~ 3 keV that cannot be accounted for by the 
Galactic absorbing column alone. This spectral model is sim¬ 
ilar to that used by Magdziarz et al. (1998) to fit these same 
data. The difference is that we do not model the so-called 
“soft excess” since we do not find significant residuals below 
~ 2.5 keV that would indicate the need for a separate spectral 
component at these energies.^ Nonetheless, we find very sim¬ 
ilar X-ray spectral indices and fluxes in the 2-10 keV band as 
do Magdziarz et al. (1998). The parameters from the fits to the 
Ginga data along with the corrected optical and UV fluxes are 
given in Table Al. 

The simultaneous EUVE, ASCA, and RXTE observations 
took place in June-August 1998 and are described in detail in 
Chiang et al. (2000). Eor uniformity’s sake, we have used only 
spectral parameters obtained from fitting RXTEfPCA data. The 
RXTE observations covered a greater period of time than did the 
ASCA observations and overlapped more with the optical obser¬ 
vations we discuss below. The PCA response matrices used in 
the original spectral analyses of these data were generated by 
the ETOOL utilities PCARMF v3.5 and PCARSP v2.36 (Chi¬ 
ang et al. 2000). However, the most recent versions of these 
programs as of September 2002, PCARMF v8.0 and PCARSP 
v8.0, produce response matrices that are sufficiently different 
from the earlier ones so that the newly fitted spectral indices 
differ from those presented in Chiang et al. (2000) by amounts 
that are ~ 3^ times larger than the 90% confidence level sta¬ 
tistical errors. In order to illustrate the effect of the systematic 
uncertainties implied by these changes, we have performed our 
analyses using the spectral parameters obtained for both sets 
of response matrices. In Table A3, the 2-10 keV fluxes and 
spectral indices are given for the “1999 PCA responses” and 
the “2002 PCA responses”. The former are taken directly from 
Table 6 of Chiang et al. 2000, while the latter have been deter¬ 
mined for this paper. 

In Table A3, the Compton reflection fractions, R, for each 
epoch are also presented. Eor the RXTE data, there is clearly 
a slight correlation between the reflection fraction and the pho¬ 
ton spectral index such that the harder power-law spectra for the 
new fits are accompanied by weaker contributions from Comp- 

which we do not use here. Only one epoch (7) of the Ginga data shows significant 
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ton reflection. However, the changes in R are not sufficiently 
large to account for the changes in the underlying spectral in¬ 
dex. 

Simultaneous UV fluxes do not exist for the 1998 
EUVEIASCAIRXTE observations, but we do have optical data 
at A5100A for seven of the epochs which were kindly provided 
to us by M. Dietrich (Dietrich et al. 2001). These data were ex¬ 
tinction and starlight corrected as described in Magdziarz et al. 
(1998). The corrected flux values are given in Table A3. Un¬ 
fortunately, the optical flux measurements were not all strictly 
simultaneous with the RXTE observations. For epoch 1.0, the 
flux listed in Table A3 was taken at JD-2450000 = 980.751, 
more than 8 hours after the RXTE observation window. For 
epochs 2.1-2.4, the optical measurements were taken either at 
the very beginning or the very end of the RXTE integration in¬ 
tervals which lasted from ^ 4 to ^ 17 hours. Given the uncer¬ 
tainties in the optical fluxes and the occasional sparseness of the 
coverage (see Dietrich et al. 2001, their Fig. 4), interpretation 
of these data as being simultaneous with the X-ray observations 
should be regarded with caution. 


3. SEMI-ANALYTIC CALCULATIONS 


As in Papers I and II, we consider a fairly simple model for 
the central regions of NGC 5548. A razor thin, flat accretion 
disk produces thermal emission that comprises the observed 
optical and UV continuum; and a centrally located, spherical 
Comptonizing region consisting of a plasma of hot electrons 
produces the power-law-like X-ray continuum detected at en¬ 
ergies ~ 0.1-100 keV. The disk is truncated at an inner radius 
I'min > lOGM/c^, the Comptonizing electrons are assumed to 
have a uniform temperature Tg, and the plasma sphere has a ra¬ 
dius Tj ~ Tniin. Thermal emission from the disk enters the sphere 
and provides seed photons for the Comptonization process. In 
turn. X-rays from the sphere illuminate the disk, heating it, and 
thus causing it to produce thermal emission in excess of that 
resulting from local viscous dissipation. The feedback between 
the disk and sphere regulates the equilibrium temperature of the 
Comptonizing electrons, and variations in the geometry through 
changes in r^in and/or determine the spectral properties of 
the X-ray, optical, and UV continua. This sort of geometry for 
thermal Comptonization in active galactic nuclei (AGNs) and 
Galactic X-ray binaries has been explored in detail by a number 
of authors (e.g., Shapiro, Lightman, & Eardley 1976; Poutanen, 
Krolik, & Ryde 1997; Dove et al. 1997; Zdziarski, Lubihski, & 
Smith 1999). 

We model the thermal emission from the surface of the ac¬ 
cretion disk as 

Fdisk — TVisc 3 " ( 1 Il)FinC5 ( 1 ) 


where a is the disk albedo. The local flux produced by viscous 
dissipation is 


3 GMM 
Stt 



( 2 ) 


where r/ is the radius of the innermost stable circular orbit. 
Note that use of this expression assumes that the torque at the 
thin disk truncation radius rmin is the same as that for a stan¬ 
dard, untruncated thin disk with zero torque at rj. We are also 
assuming that the disk is stationary, except for the variable heat¬ 
ing resulting from the incident X-rays. This incident X-ray flux 
upon the disk is given as a function of radius by 


16712 r2^(rj’ 


(3) 


where Ljc is the X-ray luminosity of the Comptonizing plasma, 
and the auxiliary function h{r/rs) is given by Zdziarski et al. 
(1999; see also Paper II). This expression assumes that the X- 
ray emission is uniform, isotropic, and optically thin throughout 
the plasma sphere. The disk emission as seen by an observer at 
a luminosity distance di and an inclination i is given by 

^ 4TTCOsihv^ 1'°° rdr 

4 c2 X„i„ /TOi[exp(/ii//kBrdisk/coi)-l]’ 
where the disk effective temperature is Tdisk = (Tdisk/crfi)*/"^, 
and /col = /coi(7disk) is a temperature dependent color correc¬ 
tion appropriate for AGN disks (Paper II; see also Hubeny et 
al. 2001). In calculating di, we assume a Hubble parameter of 
Hq = 67 km s“' Mpc“'. 

The function h{r/rs) can be approximated by 



for r > Ts (Zdziarski et al. 1999; see also Paper I), and there¬ 
fore the disk thermal emission follows an ~ r~^ distribution 
whether it results from viscous dissipation or from reprocess¬ 
ing of the X-rays from the central plasma sphere, provided that 
fs ~ I'min ^ 10 r/. Under these circumstances, the disk surface 
temperature can be modeled as /usk — To(r/■ Using 
this approximation, we have obtained estimates of the disk in¬ 
ner truncation radius rmin and the temperature To at this radius 
by fitting the optical and UV data with equation (4). Applying 
this procedure to the nine epochs of the 1989-90 data, we have 
found the rmin-values listed in Table A2. These results depend 
sensitively on the disk inclination, and that dependence affects 
the inferred values for the other parameters describing the X- 
ray portion of the SEDs. We will explore the implications of 
this dependence in §5. 

The underlying relationship between the X-ray luminosity, 
the black hole mass, and the accretion rate can be seen by 
using the approximation for h(r/rs) (eq. [5]) and the afore¬ 
mentioned fits for r^iin and Tq. Since broad line region rever¬ 
beration mapping estimates of the central black hole mass in 
NGC 5548 give M,™ = 6 x lO’ M© (Peterson & Wandel 2000), 
and r^iin > 20 for all nine of the 1989-90 epochs, we 

can neglect the term containing rj jr in the expression for Tvisc 
(eq. [2]). Using these approximations, the black hole mass and 
the accretion rate are related to the X-ray luminosity and the 
radius of the plasma sphere by 

GMM= (6) 

8 3 

Since we are assuming that the disk is stationary (eq. [2]), the 
non-varying part of the opticaFUV SED is modeled by the 
MM term on the left-hand side of equation (6), while any op- 
tical/UV variability in the model must be caused by changes 
in the thermal reprocessing component represented by the 
term and/or by changes in the truncation radius rmm- 

The values of and are also constrained by the observed 
2-10 keV X-ray flux and spectral index, and by two energy 
balance conditions. The first energy balance condition follows 
from an empirical relation between the X-ray spectral index and 
the Compton amplification factor A, which is the ratio of the 
X-ray luminosity and the luminosity of “seed” photons, Lj, en¬ 
tering the Comptonizing plasma: 

r = ro(A-i)-'/‘5 (7) 


:r„|U-, 


-1/s 


( 8 ) 
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The parameters To = 2.15 and S = 14 have been found from fits 
to Monte Carlo calculations of the thermal Comptonization pro¬ 
cess (Malzac, Beloborodov, & Poutanen 2001; see also §4). If 
the seed photons consist only of thermal emission from the disk, 
then Ls is determined by the disk flux distribution Fdisk(?‘) and 
by the geometry of the disk and Comptonizing plasma through 
rs and rmin • 

pOO 

Ls = 4Tr / g(r,f,)I(r)rdr, (9) 

J rmin 


where/(r) = Fdisk(r)/7r is the disk surface brightness, f = 
and g(r,rs) is defined in Paper II and equals the fraction of 
isotropic emission that enters the central plasma sphere from 
a disk annulus (f,r + dr). Modeling the X-ray spectrum as an 
exponentially cut-off power-law, Le cx E'^~^ exp(-£’ /£e), where 
Ec = 2kBTe, a value for can be found using the 2-10 keV X- 
ray continuum flux, F 210 , and X-ray spectral index, F: 


Lx = ^4TrdfF2io 


_ Q(r,E^ijE,) _ 

e(r,2keV/£,)-e(r, 10keV/£,)' 


( 10 ) 


Here Emin = 0.01 keV is an assumed nominal lower cut-off en¬ 
ergy, the precise value of which is not important for photon 
indices F < 2, and 


Q(T,x)= / (11) 

J X 

Note that this expression is the incomplete gamma function. 
The factor ^ in equation (10) parameterizes any anisotropy in 
the thermal Compton emission. If this emission were isotropic 
and the plasma sphere were optically thin, then $ = 1. Be¬ 
cause photons from the first scattering order are preferentially 
scattered back towards the disk, the X-ray intensity incident 
upon the disk can be larger than the intensity seen by a dis¬ 
tant observer (Paper II). Our Monte Carlo calculations show 
that ^ ~ 1.5 does a reasonable job of modeling this anisotropy 
for observer inclinations i ~ 30^5° (see §4). 

The second energy balance condition is more global in na¬ 
ture and comes from the expectation that the hot plasma is less 
efficient at radiating accretion power than a standard optically 
thick accretion disk would be in its place: 


Ly+ Ltj 


( 12 ) 


The quantities Ltr and Luntr are the luminosities resulting from 
viscous dissipation for a disk truncated at an inner radius rmin 
and an “untruncated” disk extending to the innermost stable cir¬ 
cular orbit r/, so that 

pOC) 

i-muntr = 47r / r(fr/Visc(r). (13) 

J rm\n,Ll 

For a given value of r/, equation (12) is equivalent to 

Lx<7jMc^-L^, (14) 

where the efficiency p = 0.5 for r/ = GM/c^f' We generally as¬ 
sume that a firm upper limit to the X-ray luminosity is given by 
?7 = 1. Of course, this high a value would correspond to accre¬ 
tion efficiencies that exceed both that of the standard thin disk 
and those expected for other models such as advection domi¬ 
nated accretion flows (ADAFs; Narayan & Yi 1995). However, 
even higher efficiencies are conceivable if the spin energy of 
the central black hole can be tapped (e.g., Krolik 1999; Gam- 
mie 1999; Wilms et al. 2001). 


To summarize, we effectively have five relations constraining 
the properties of the inner accretion flow. We list them here in 
terms of their functional dependence on the model parameters 
and the observationally determined quantities: 


MM 

— AlAl (Lx min •> ) 5 

(15a) 

F 

= r(L^,L^;ro,(5), 

(15b) 

Ls 

— Ls(Lx’)f's-i^mm')^^^') 

(15c) 

Lx 

= L,(r,;F2io,r), 

(15d) 

Lx 


(15e) 


The first relation represents equations (1^) (or approximately 
eq. [6] by itself); and the remaining four relations represent 
equations (8), (9), (10), and (14), respectively. In these five 
relations, there are six unknowns^. Lx, Ls, Tg, rs, MM, and ri', 
and we consider the remaining quantities as being either di¬ 
rectly measured from or constrained by the observations: F 210 
and F from the X-ray data, and r^ain by the optical and UV data 
as represented by Fbuv (see eq. [4]). In order to obtain a unique 
solution to equations (15), additional information is necessary. 
Direct measurement of the roll-over energy Ec would provide 
the most useful constraints since that would provide estimates 
of both Lx and Tg, but instruments with the required sensitiv¬ 
ity in the relevant energy range are not yet available. Rever¬ 
beration mapping measurements of the black hole mass can be 
used to estimate r/; but as we will see below, the standard value 
for NGC 5548, = 6 x 10’ Mq (Peterson & Wandel 2000), 

yields accretion rates so low that equation (14) is violated or 
that imply electron temperatures that are inconsistent with the 
observed X-ray continua. 


4. RESULTS 

In order to get a sense for how the quantities rj and MM can 
affect the implied solutions, we have solved equations (15a-d) 
numerically over a grid of values spanning M = (1-7) x IO’Mq 
and MM = (0.5-7.5) x 10^ Mq’ yr“\ assuming rj = GM/c^. 
As an example, in Figure Al, the values of ksTe obtained from 
these fits are plotted versus MM for the 1989-90 epoch 2 data. 
Since rmin = 1.8 x 10'“* cm > lOr^ over the entire range of 
masses, the r//r term in equation (1) plays a minor role in de¬ 
termining the disk flux distribution, and equation (6) roughly 
applies. Consequently, the electron temperatures are nearly in¬ 
dependent of M and are determined mostly by the quantity MM. 
This functional dependence is shown in Figure Al by the var¬ 
ious symbols (corresponding to the different masses) and by 
the single solid curve that connects the values of ksTg aver¬ 
aged over the masses at each MM value. For each black hole 
mass, we compute the minimum accretion rate that satisfies 
equation (14), 

. V Lx + Lxx 

Mpcim = J , (lb) 

r]c^ 

where we have taken the accretion efficiency to be p = 0.5. The 
dashed curves show the fitted kgTc values plotted versus MWmin 
for the specified mass; solutions of equation 15a-d that also 
obey equation 15e will lie to the right of the dashed curves. 
Hence, the intersection of each dashed curve with the solid 
curve gives the minimum value of MM that satisfies the en¬ 
ergy balance condition (eq. [14]). That intersection point also 


This is obtained from equation (13). We do not use the more familiar value of r; = 0.42 since we are neglecting general relativistic effects throughout this work. 

^ In fact, the disk inclination i is another unknown. For clarity, we defer explicit discussion of it until §5 and assume cos / = 0.7 for our analyses of the NGC 5548 data 
unless otherwise indicated. 
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corresponds to the maximum electron temperature that can be 
achieved for the given mass. ForM = 10^ M©, the electron tem¬ 
perature during epoch 2 must satisfy hgTe < 200 keV; while for 
M = 6 X 10 ^ Mq, it must satisfy ksTg < 60 keV. 

Several observations of the same source in different flux 
states can place additional limits on the black hole mass M and 
accretion rate M, provided the latter is fixed as we are assuming 
by adopting a stationary thin disk structure. As we described 
above, we use the optical and UV continuum fluxes to deter¬ 
mine the disk inner radius, rniin, in each epoch. In contrast to 
our previous studies of NGC 3516 and NGC 7469 (Paper II), 
the optical and UV data for NGC 5548 do seem to require vari¬ 
ations in Tniin. This is evidenced by the relatively large opti- 
cal/UV color changes occurring over these observations (Ta¬ 
ble Al). These observations are separated in time by intervals 
ranging from two days to nearly one year. If variations in r^in 
result from viscous processes, then any feasible observational 
time scale would be far too short to expect significant changes. 
However, if the rmin-variations occur because of thermal insta¬ 
bilities causing the disk to alternately “puff up” and then “de¬ 
flate” near the innermost radius of the thin disk, then such vari¬ 
ations could occur on a thermal time scale: 


fth ~ 10 S 


(f^) 


-1 


.3/2 


M 


- 1/2 


(17) 


,2xl0'4cmy VIO^Mq, 
where a is the anomalous viscosity parameter. This is roughly 
consistent with the time scales of the rmin-variations listed in Ta¬ 
ble A2, except perhaps for the 2-day interval between epochs 4 
and 5 that would require somewhat larger values of a and/or 
black hole mass than are indicated in equation (17). 

In Figure A2, we plot the k^Tg-MM curves for all nine 
of the 1989-90 epochs. The vertical dotted lines indicate 
the minimum MM values for epoch 1 assuming masses M = 
(1,2,3, 4, 5, 6, 7) X 10^ M©, ordered left-to-right. These lower 
limits are greater than the corresponding ones found for the 
other eight epochs, and so they represent the most restrictive 
minimum values for MM. Given these lower limits, a black 
hole mass of /W = 6 x 10^ Mq implies a minimum value of 
(MM)niin — 7.5 X 10^ Mgyr"^ and requires rather low electron 
temperatures for the eight other epochs, with ksTg < 20 keV 
in all cases. For temperatures that low, the thermal roll-over 
would be detected in the X-ray spectra taken by instruments 
such as Ginga and RXTEIVCA, which have useful sensitivity 
up to ~ 20-30 keV for the brighter type 1 Seyfert galaxies. Of 
course, the degree to which these observations could detect the 
effects of a low electron temperature depends on the shape of 
the thermal Comptonization continuum. Although power-laws 
with exponential cut-offs are often used to model these spec¬ 
tra, detailed thermal Comptonization calculations, such as our 
Monte Carlo calculations (§ 5), generally show sharper high en¬ 
ergy cut-offs so that the thermal roll-over may not be as evident 
at lower energies. Furthermore, a strong Compton reflection 
component can hide a low energy thermal roll-over much in 
the same way it can mask a softer underlying power-law in¬ 
dex. Nonetheless, one can adopt a fairly conservative mini¬ 
mum value for 7), as an observational constraint. For our pur¬ 
poses, we will require ksTg > 30 keV, which is indicated by the 
lower dotted horizontal line in Figure A2. Since the ksTg-MM 
curves decline fairly rapidly, a much lower minimum temper¬ 
ature will not drastically alter the implied limits on M and M. 
From Figure A2, we see that the epoch 7 and 9 curves con¬ 
strain MM < 2.5 X 10^ Mgyr"' and the black hole mass to be 
M < 2 X 10^ Mq. Since we expect the accretion efficiency 


to be less than maximal, we will take M = 1 x 10^ Mq and 
M = 2 X 10“^ Mq yr“^ as nominal values hereafter. Assum¬ 
ing r/ = GM/c^, this yields Luntr = Mc^/2 = 6 x 10"^“* erg s“^ as 
our heuristic upper-limit for the bolometric luminosity of this 
source (eq. [14]). Note that this is 60 percent of the Edding¬ 
ton luminosity at the chosen black hole mass, which would be 
difficult to reconcile with an ADAF interpretation for the inner 
accretion flow. 

The model parameters and SEDs which we find for the 1989- 
90 data using these values are listed in Table A2 and shown in 
Eigure A3, respectively. In all cases we find that rj/rn,in ^ 1, 
with this ratio being slightly greater or slightly less than unity 
at different epochs. Note that values less than unity need not 
imply that there is a physical gap between the hot plasma and 
the inner truncation radius of the disk. Instead, it may indicate 
a central accretion flow that is more flattened and aspherical 
compared to the sphere we assume in the present model (see 
Paper II). 

The SEDs shown in Eigure A3 often have more power in 
the Comptonized X-ray component than in the opticaFUV disk 
thermal component. This is still consistent with a radiatively in¬ 
efficient hot inner flow, because that efficiency should be com¬ 
pared to a radiatively efficient disk at the same radius where 
most of the accretion power is generated. Indeed, the very fact 
that the X-rays appear to take up such a large fraction of the 
bolometric power is strongly suggestive of the accretion flow 
geometry we are considering here. 

Applying our calculations to the 1998 observations is prob¬ 
lematic since we lack the UV data required to constrain the disk 
inner radius. However, having ascertained M and M from the 
1989-90 data, we can explore the feasible range of disk inner 
radii for the 1998 epochs in light of our temperature constraint. 
In Eigure A4, we plot the equilibrium electron temperature ver¬ 
sus Tmin for the seven epochs of the 1998 data for both sets of 
response matrices. For clarity, we have plotted the individual 
curves that are obtained from analyses using the 1999 PCA re¬ 
sponse matrices; while for the 2002 response matrices, we sim¬ 
ply show a hashed region that is bound by the minimum and 
maximum ksTg values at each value of rn,in. 

Using the 1999 response matrices, the curve for epoch 2.2 
requires that r^ain < 2.4 x 10^"^ cm (= 160 G(M/10 ^Mq)/c^) in 
order to satisfy ksTg > 30 keV. If we assume that this upper 
limit for rn,in applies to all of the 1998 observations, then the 
electron temperature for epoch 3.2 would be ksTg > 1500 keV 
(see Table A4). Since we do not consider pair-production or the 
possibility of a non-thermal component of relativistic electrons, 
another constraint we are tacitly imposing in this model is that 
the electron temperatures should satisfy kgTe < nigC^. For tem¬ 
peratures much greater than this, the thermalization time scales 
exceed the cooling time scales, and a non-thermal electron com¬ 
ponent will arise (Coppi 1999). The presence of a signifi¬ 
cant non-thermal component would invalidate the energy bal¬ 
ance implied by the thermostatic interaction between the disk 
reprocessed radiation and the Comptonization emission from 
the hot plasma. Therefore, in order to satisfy this implicit up¬ 
per limit for the temperature, the curves shown in Figure A4 
suggest that the disk inner radii for these data span a range 
at least as wide as r^in — (2.4-3.2) x 10^"^ cm. However, for 
the spectral parameters derived using the 2002 response matri¬ 
ces, the disk inner radius can take any single value in the range 
(3.1-4.7) X 10'“* cm for all seven epochs and still be consistent 
with electron temperatures that lie within our nominal range of 
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30-511 keV. Since ~ for all these calculations (see Ta¬ 
bles A4 and A5), it is worth noting that, despite the changes in 
the RXTE/PCA response matrices, all of these radii are roughly 
consistent with our estimates of the size of the Comptonizing 
region that we derived from the temporal lags between the soft 
(0.16 keV) and hard (2-10 keV) X-ray light curves of the 1998 
EUVE/ASCA/RXTE observations (Chiang et al. 2000). 

In Figure A5, we plot model spectral energy distributions 
for the seven 1998 ground-based/T^XTiE epochs. Two of the 
sets of curves have been computed using the parameters from 
the 1999 response matrix spectral fits, assuming r^in = 2.4 x 
10^"^ cm (solid curves) and rmin = 3.2 x 10'“* cm (dashed). The 
SEDs corresponding to the 2002 response matrix fits are also 
shown (dot-dashed curves) and have been computed assuming 
^■min = 3.2 X 10'“* cm. The disk luminosities resulting from vis¬ 
cous dissipation are Ltr = 1.0 x 10"^^ and 0.7 x 10"^^ erg s“^ for 
^■min = 2.4 X 10^"^ cm and 3.2 x 10^"^ cm, respectively. Global 
energy balance, equation (14), is obeyed in all cases. The SEDs 
shown in Eigure A5 illustrate the dramatic effect that the change 
in the response matrices have on the inferred broad band X- 
ray continuum. It is worth emphasizing, however, that the 2- 
10 keV continuum fluxes obtained using the two different sets 
of response matrices are fairly consistent with one another (see 
Table A3). Consequently, the largest differences in the model 
SEDs occur at energies that lie well outside the range directly 
measured by RXTE/PCA. 

5. MONTE CARLO CALCULATIONS 

We have implemented a three dimensional Monte Carlo cal¬ 
culation of the thermal Comptonization/thermal reprocessing 
model we have described in the previous section. A more de¬ 
tailed description appears in the Appendix and in Chiang & 
Blaes (2002). In this section, we compare results from our 
Monte Carlo calculations with the empirical relations we have 
used in the previous section and with the spectral energy distri¬ 
butions obtained using those semi-analytic methods. 

In Eigure A6, we plot the X-ray continuum spectral index 
versus the amplification factor, A = Ly^jLs. The individual data 
points are from our Monte Carlo calculations, the solid curve is 
the parameterization of Malzac et al. (2001; see eq. [8]), and the 
dashed curve is a best fit to the Monte Carlo data with Eq = 2.06 
and 8 = 15.6. The filled triangles represent the Monte Carlo cal¬ 
culations that produced the histogrammed SEDs shown in Eig¬ 
ure A3; we discuss these simulations in greater detail below. 
The other symbols represent simulations that were computed 
for a black hole mass of M = 2 x IO^Mq, a disk inner radius 
^min = 2 X 10'“* cm, an observer inclination of / = 30°, five mass 
accretion rates M = 0,0.01,0.02,0.04, and 0.08 Mq yr“', and 
five sphere radii = 0.5,0.75,1.0,1.25, and 1.5 x 10'“* cm for 
a total of 25 spectra. 

In Eigure A7, we plot the ratio Ej,repr/Tj:,app versus 
for these same two sets of simulations. Here Lj^repr is the seed 
photon luminosity produced by thermal reprocessing (i.e., the 
total Es minus the portion produced by viscous dissipation), 
and Ex.app is the apparent thermal Comptonization luminosity, 
i.e., the luminosity one would compute from equation (10) for 
^ = 1. The solid curve is the relation one would find if the ther¬ 
mal Comptonization luminosity were indeed isotropic so that 
the flux incident upon the disk is given by equation (3). The 
dashed curve is the solid curve multiplied by ^ = 1.5, the value 
we adopted to model the anisotropy of the thermal Comptoniza¬ 
tion emission. This factor also takes into account differences 


in the albedo of the disk. Eor the semi-analytic calculations, 
we assumed a constant albedo of a = 0.15, while the disk albe¬ 
dos found in the Monte Carlo calculations vary with the inci¬ 
dent spectrum, ranging from a ~ 0.2 for E ~ 1.9 to a ~ 0.4 for 
E ~ 1.5. As can be seen from Figure Al, for r^/rnim ^ 1.1, a 
value of ^ = 1.5 is appropriate, while for lower values, which 
correspond to weaker coupling between the hot plasma and 
cold disk and therefore produce harder spectra, a value of ^ ~ 1 
should probably be used. 

In Figure A3, the histograms are the Monte Carlo SEDs we 
computed using the values of M, M, r^, and r^in found in §3 for 
the 1989-90 data. The equilibrium electron temperatures we 
obtain from our Monte Carlo calculations are systematically 
larger by ~ 20% than those from our semi-analytic methods. 
Now, in the Monte Carlo simulations, we have set the radial 
Thomson depth, r,, of the Comptonizing sphere to be equal to 
the Thomson depth, r, found from the observed spectral index, 
the electron temperature from the semi-analytic method, and 
the empirical relation (Beloborodov 1999), 

r = (18) 

where the Compton y-parameter is 

y = T(l+r)40,(l+40,), (19) 

and 9e = kETelnigC^. Since the spectral index is set by the 
geometric factors, and r^in, E is independent of the Thom¬ 
son depth of the plasma. However, the amplitude of the ther¬ 
mal Comptonization continuum is determined by the amount 
of seed photons scattered in the Comptonizing region and so 
should be roughly proportional to 1 -e~'^'. Therefore, the mis¬ 
match between the continuum levels of the semi-analytic and 
Monte Carlo calculations, most clearly seen in the epoch 5 re¬ 
sults, can be reduced by adjusting the radial Thomson depth 
of the plasma sphere. This will, of course, result in different 
electron temperatures as well. 

For the epoch 1 data, the shape of the Monte Carlo SED is 
more difficult to reconcile with the exponentially cut-off power- 
law used in our semi-analytic method because for sufficiently 
high electron temperatures the individual scattering orders are 
partially resolved in the X-ray spectra. This illustrates another, 
observationally-based motivation for an upper limit on the elec¬ 
tron temperature. If the mean separation in energy between the 
peaks of individual scattering orders exceeds their width, these 
orders will be partially resolved in the resulting thermal Comp¬ 
tonization continuum and will appear as a number of spectral 
breaks and/or local maxima and minima (e.g., Skibo & Dermer 
1995). Empirically, we find from our Monte Carlo calculations 
that an upper limit of < nigC^ is required to prevent these 
features from appearing. 

6. DISCUSSION 

The wide range of electron temperatures, ~ 40-700 keV, 

that we found for the ground-based/Zt/E/G/nga data is driven by 
the relatively broad range of opticaEUV fluxes (see Table Al). 
In Paper II, we obtained similar results for the NGC 7469 moni¬ 
toring observations, which also had large optical/UV variations. 
In that work, we were able to narrow the range of temperature 
changes that was required by allowing the X-ray spectral in¬ 
dices to vary within their 1-a uncertainties. Eor the NGC 5548 
ground-based//t/E/G/ngfl data, the most conspicuous out-lier is 
the fit to the epoch 1 data. Those data yield an electron tem¬ 
perature that is at least three times larger than the temperatures 
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found for the other eight epochs. The epoch 1 data also have 
an unusually hard X-ray spectral index, T = 1.56. This value 
lies at the very low end of the range of spectral indices that 
are typical for type 1 Seyfert galaxies (Zdziarski, Poutanen, & 
Johnson 2000). Furthermore, it is substantially smaller than the 
indices of the other Ginga epochs (with the notable exception of 
epoch 8). However, as we noted in Paper II, softer spectra have 
less X-ray flux at energies >10 keV than do harder spectra for 
a given value of hgTe and a fixed 2-10 keV flux. In order to re¬ 
duce kgTg for epoch 1, we would either have to make the X-ray 
spectrum even harder, i.e., use F < 1.56, or take a larger value 
of the 2-10 keV flux. Setting T 210 = 3.73 x 10“" erg cm“^s“' 
for epoch 1, which corresponds to a 10% increase in 2-10 keV 
flux over the fitted value (see Table Al), we obtain a uniform 
20% reduction in the equilibrium electron temperature for all 
values of MM considered. The resulting changes in the inferred 
minimum accretion rates are even smaller: forM = 7x10^ M©, 
Mmin is reduced by 13%, while for M = 2x 10^ M©, the reduc¬ 
tion is less than 4%. 

For the preceding analyses of the NGC 5548 data we have 
assumed an observer inclination of i = 45°.5 (cos / = 0.7). Our 
calculations are surprisingly sensitive to this value. Since the 
optical/UV spectral shape for a given epoch is determined 
largely by the temperature of the disk at the innermost radius, 
7b is independent of cos/. The disk flux is proportional to 
~ ^ givcn set of opticaFUV data, we 

have Tmin oc (cos/)“'/^. If the disk emission is dominated by 
the thermally reprocessed component, then the geometry is de¬ 
termined by the X-ray spectral index (Zdziarski et al. 1999), 
and we have oc oc (cos/)“'/^, Lj oc oc (cos/)“', and 

Lx oc (cos/)“'. Inserting these relations into equation (6), we 
obtain 


87raji7bV^i„Q L^,or,^o(l-a) 
3(cos/)^/2 8(cos/)^/^ ’ 


( 20 ) 


where Lx.q, Tmin.o, and r^o are fiducial values evaluated for 
cos / = 1. The implied values of MM will therefore be larger 
for smaller cos /. Similarly, the global energy balance condition 
(eq. [14]) allows us to obtain a (cos/)-scaling forM: 


4wB7oV^i 0 Lx,q 

rjMc > -^-l-[l-(l-a)/(]— 

cos/ cos/ 


( 21 ) 


where R is the fraction of the X-ray luminosity intercepted by 
the disk. Taking the minimum value for M in order to have a 
maximum value for the black hole mass, this expression and 
equation (20) imply M oc (cos/)“'/^ and M oc (cosi)~K Now, 
the X-ray observations and our ad hoc limits for the electron 
temperature fix the range of thermal Comptonization roll-over 
energies, and therefore they essentially fix the Lx values for the 
1989-90 data as a whole (see Figs. A2 and A8), so the second 
term on the right-hand side of equation (20) will, in practice, go 
as (cos/)“'/^. Nevertheless, the general trend forM andM to be 
larger for smaller values of cos / still holds. In fact, we find for 
these data that each term on the right-hand side of equation (20) 
is much larger than the left-hand side, so given the different ef¬ 
fective scalings with cos / of those two terms, the value of MM 
is quite sensitive to the assumed inclination. 

Figure A8 illustrates more concretely how much the feasible 
ranges of M and M depend on the choice of the observer incli¬ 
nation; we show plots of ksT^ versus MM for / = 30° and / = 60° 
for the 1989-90 data. If we take our minimum electron temper¬ 
ature limit seriously, then for / = 30°, we find M < 2 x 10® M© 
and MM < 2 x 10“* M^yr"'. By contrast, for / = 60°, we find 


M < 3 X IO^Mq and MM < 6 x 10® M^yr'^ In the latter 
case, an electron temperature of kgTe ~ 1200 keV is implied 
for epoch 1, which well exceeds our upper temperature limit. 
For the former case, although the electron temperatures fall 
within our adopted limits of 30 keV < ksTe < nigC^, the implied 
black hole mass is a factor of 30 smaller than the reverbera¬ 
tion mapping estimate. Therefore, as a compromise, we have 
taken / ~ 45° as our nominal value which does yield a some¬ 
what high electron temperature (ksTg ~ 700 keV) for epoch 1, 
but also gives a mass upper limit (M < 2 x 10® Mq) that is not 
in too great conflict with the reverberation mapping estimates 
(see Krolik 2001). Note that extraction of black hole spin en¬ 
ergy, which we have not considered here, would ease the upper 
limits on the black hole mass by allowing, in principle, accre¬ 
tion efficiencies in excess of unity. 

Even assuming a mass as large as the reverberation map¬ 
ping estimates, our fits to the optical/UV continua give r^in > 
20GM/c®. This is marginally inconsistent with the range of 
disk inner radii of rmin = 7.5-15 GM/c® found by fitting a rela- 
tivistically broadened iron line to the 1998 ASCA data (Chiang 
et al. 2000; see also Mushotzky et al. 1995). IfM < 2 x 10® M©, 
then our results would be in serious disagreement with the de¬ 
tection of a relativistically broadened line. This conclusion re¬ 
lies on the 1989-90 epoch 1 data providing the relevant limits 
on MM. It is clear from Figures A2 and A3 that the SED dur¬ 
ing this epoch was unusual, and the source may have been in a 
transient state for which the X-ray luminosity briefly exceeded 
the limit imposed by equation (14). Averaging over the nine 
1989-90 epochs, we obtain a much lower set of MM-limits than 
we find from epoch 1 alone. These are indicated in Eigure A2 
by the arrows along the top of the plot. We still require that 
MM < 2.5 X 10® Mgyr”' in order to have reasonable electron 
temperatures for epochs 7 and 9; and so these averaged limits 
imply M < 3 X 10® Mq. If we take M = M^ = 6 x 10® Mq, 
then M < 4.2 x 10“® Mq yr“' and Mc®/2 = 1.2 x 10"^^ erg s“'. 
In this case, the X-ray luminosity exceeds the limit imposed by 
global energy balance for all observations except epoch 9. 

If the Comptonizing plasma is an ADAE, then the radial ve¬ 
locity of the accretion flow in this region is Vr ~ av^, substan¬ 
tially higher than that of the surrounding thin disk (Narayan & 
Yi 1995). Erom mass conservation and our estimates of the ac¬ 
cretion rate, we can determine a characteristic Thomson depth, 
Tc, for the central region. Mass conservation in an accretion 
disk is given by 

M = 27rrVrS (22) 

~ 47rr®pvr, (23) 


where S = 2hp is the disk surface density, h is the half¬ 
thickness, and we have set h ^ r for the central region which 
we assume to be dominated by ions at the virial temperature. 
The density of the central region is roughly 


/3~ 


Tc 

fsKe" 


(24) 


where Ke = 0.4 cm^g“^ is the electron scattering opacity. Solv¬ 
ing equations (23) and (24) for the Thomson depth and inserting 
nominal values, we obtain 


Tc ~ 0 


-iit) 


-1 


- 1/2 


2 X 10^"^cm^ 

M mm \ 

2x 10 ®Mq ) V23^^T0®M|y ^) ■ 



Despite the crudeness of this estimate, it is encouraging that Tc 
lies within the range of Thomson depths yielded by our model 
(Tables A2, A4, and A5). 

A potentially serious problem for this model is the large 
change in Thomson depth that we find between epochs 2.1 
and 2.2 of the 1998 observations for both sets of response ma¬ 
trices (Tables A4 and A5). Since changes in r^ain must occur 
on time scales that are at least as long as the thermal time scale 
and since epochs 2.1-2.4 are fairly closely spaced over a 240 ks 
period (Chiang et al. 2000; see also Table A3), we must adopt a 
single radius for these epochs. For the 1999 response matrices 
and assuming r^in = 2.4 x 10'“* cm, the Thomson depth must 
change from t = 0.06 to t = 2.26 over ~ 30 ks. Even for the 
2002 response matrices, the Thomson depth must change from 
T = 0.21 to 0.69, assuming rmin = 3.2 x 10^"^ cm. The necessity 
for these large changes can be readily understood. During these 
epochs, the A5100A flux was essentially constant (Table A3). 
As we found for NGC 3516, in order to have a nearly con¬ 
stant optical flux, the radius of the Comptonizing region has 
to maintain an inverse relationship with the X-ray luminosity 
(eq. [6]). Since the X-ray spectral index increased over this in¬ 
terval, the radius also had to increase in order to produce a 
softer spectrum; and therefore had to decrease. However, 
the 2-10 keV X-ray flux increased at the same time, and so the 
electron temperature had to decrease by a very large amount 
to compensate. For a fixed Thomson depth, the electron tem¬ 
perature after the transition would result in an X-ray spectrum 
that would be much softer than is observed, and so the Thom¬ 
son depth had to change substantially as well. The time scale 
for the evolution of the Thomson depth of the central region 
should be limited by either the thermal or accretion (i.e., vis¬ 
cous) time scales. Assuming a virial temperature for the ions in 
the central region, these two time scales are approximately the 
same, 4cc ~ ~ 10® s (eq. [17]). 

One possible mitigating factor here is the relatively sparse 
coverage of the optical observations (see §2). Multiplying the 
optical fluxes for epochs 2.1-2.4 in Table A3 by the factors 
0.86, 1.14, 1.07, and 0.98, respectively, yields a nearly constant 
Thomson depth of t = 0.44 and electron temperatures hgTe ~ 
169-184 keV, assuming rmin = 2.4 x 10'“* cm. Unfortunately, a 
> 28% change in the A5100A flux over the ~ 30 ks separat¬ 
ing epochs 2.1 and 2.2 does not seem to be supported by the 
data (Dietrich et al. 2001, see their Fig. 4). Of course, it is the 
UV photons produced near the innermost radius of the thin disk 
that are most important for determining the energy balance of 
the Comptonizing plasma, and they may not be well correlated 
with the emission at longer wavelengths on these time scales. 
The disk flux at A5 lOOA is produced at a characteristic radius of 
Oioo ~ rn,in[7’o/(0.29 cm K/5.1 x 10“® cm)]^/^ ^ 3 x 10^® cm. 
Therefore, the optical emission should be smeared out and de¬ 
layed relative to the UV emission by ~ 10® s. 

7. CONCLUSIONS 

We have fit simultaneous optical, UV, and X-ray data of NGC 
5548 with reasonable success over a wide range of continuum 
flux levels using just a crude model consisting of a stationary, 
truncated accretion disk and a spherical hot inner plasma. Un¬ 
der fairly general assumptions concerning accretion efficiency, 
energy balance, and self-consistency, we are able to place up¬ 
per limits on the black hole mass, M <2 x 10^ Mq and the 
accretion rate, M < 2.5 x 10® MQyr“'/M. The upper limit on 
the black hole mass is less than that derived from reverbera¬ 


tion mapping, which may itself be subject to some systematic 
uncertainty. Increasing the accretion efficiency, perhaps by ex¬ 
tracting black hole spin energy, would relax our constraint on 
the black hole mass. We also find that the disk truncation radius 
is of order the size of the Comptonizing plasma and is ~ 2- 
5 X 10'“* cm. This is remarkably consistent with the size scale 
inferred in a completely different manner from the measured 
EUV to X-ray time lags (Chiang et al. 2000). 

Although the present model does not provide a perfect de¬ 
scription of the multi wavelength variability of NGC 5548, these 
calculations do show the potential for using multiwavelength 
data for constraining properties of the disk/black hole system 
in these objects. Wider spectral coverage with more sensi¬ 
tive X-ray observatories will help enormously in testing models 
such as this. In particular, the ability to measure the thermal 
Comptonization roll-over at energies > 50 keV will allow us 
to abandon our ad hoc temperature limits and the assumption 
of a global energy balance condition. Forthcoming INTEGRAL 
observations are the best near-term prospects for better measur¬ 
ing the thermal Comptonization roll-over in Seyfert 1 galaxies; 
and some of the stronger sources — NGC 4151, IC 4329A, and 
MCG-6-30-15 — are among its highest priority targets. How¬ 
ever, even INTEGRAL does not have the sensitivity to measure 
changes in the roll-over energy on the required time scales of 
~ 30 ks. For this, an instrument such as the hard X-ray detector 
aboard ASTRO-E2 is required. 

XMM-Newton observations of NGC 5548 alone may be able 
to provide stringent tests of this model. If the iron Ka line is rel- 
ativistically broadened, then the EPIC detectors aboard XMM- 
Newton should be able to determine both the disk inclination i 
and the disk inner radius, rmin, in units of GM/c^. Combined 
with optical/UV spectroscopy available from the XMM-Newton 
Optical Monitor that can determine r^in directly, such observa¬ 
tions should be able to constrain the black hole mass. However, 
if I'min 3> GMjc^, as suggested by our analysis of the NGC 5548 
data, then the iron line will not be relativistically broadened, 
and XMM-Newton observations will provide at best a weak up¬ 
per limit on M. 

We wish to emphasize again that our model assumes that the 
disk itself is stationary, and all optical/UV variability is a re¬ 
sult of the variable incident X-ray flux and varying r^in. As we 
noted in equation (17), the disk can be out of thermal equilib¬ 
rium in the innermost regions on time scales comparable to the 
observed variations. If thermal instabilities are responsible for 
the variations in rmin, then our model may not be fully consis¬ 
tent. Even in the optically emitting regions, the thermal time 
scale is only ~ 2 yr (a/0.1)“' for a black hole mass of 10^ Mq. 
It may therefore be interesting to develop our model further to 
incorporate a non-stationary disk structure. 
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APPENDIX 

DESCRIPTION OF THE THERMAL COMPTONIZATION MONTE CARLO 

Seed photons from the disk are drawn assuming that the disk radial temperature profile is set by contributions from local viscous 
dissipation and reprocessing of X-ray flux from the hot plasma. The intensity of the emission is assumed to be isotropic at each point 
on the surface of the disk, and the energy distribution of the photons is modified from black body to account for the temperature 
dependent empirical color correction we have applied (see §3 and Paper II). Eor disk photons which intercept the central Comptoniz- 
ing region, the optical depth is computed using the full Klein-Nishina cross-section and assuming a single-temperature relativistic 
Maxwellian distribution for the electrons (Poutanen & Svensson 1996). A scattering location along the photon path is drawn from the 
appropriate distribution; and if the photon is scattered, then its subsequent trajectory and energy are drawn from the scalar Compton 
redistribution function (Poutanen & Svensson 1996; Jones 1968). The optical depth to the boundary of the Comptonizing region is 
computed and a new scattering location is drawn. This continues until the photon escapes the Comptonizing region. 

If the escaping photon intercepts the disk, then a separate Compton reflection Monte Carlo is employed. In this calculation, the 
incident photon strikes a semi-infinite slab of material at the appropriate angle and is scattered by the cold electrons in the slab until 
it either re-emerges or is absorbed. Eor the material composing the slab, we adopt the absorption cross-sections and abundances 
of Balucinska-Church & McCammon (1992), and assume that the hydrogen and helium atoms are fully ionized and the remaining 
elements are neutral. This part of the code is functionally identical to the Monte Carlo calculations used by Magdziarz & Zdziarski 
(1995) to derive parameterizations of the Compton reflection Green functions. We have tested our code against their expressions and 
find consistent reflection spectra. Photons that escape the Comptonizing region without hitting the disk and those that emerge from 
the disk without re-entering the plasma are included in the final spectra according to their energy and direction. 

The thermal Comptonization part of our Monte Carlo code has been checked against a two-dimensional iterative scattering (IS) 
code in which the integrals for the source functions at each scattering order are computed numerically. This IS code has in turn been 
checked against the one-dimensional IS code of Poutanen & Svensson (1996). We confirm that our Monte Carlo calculation agrees 
at each scattering order with the IS code. Using the Monte Carlo code rather than the IS code presents several advantages. The 
principal one is that it runs more quickly than the IS code, within a certain accuracy, allowing us to explore parameter space more 
efficiently. Eurthermore, interfacing it with the Compton reflection part of the code is easier since both methods deal with individual 
photons. Lastly, should we decide to generalize this calculation, it would be a straight-forward matter to modify the geometry of the 
Comptonizing region or take into account curved spacetimes in computing the trajectories of the photons. 

In order to find the equilibrium temperature of the Comptonizing electrons, we proceeded iteratively, using the difference between 
the model and observed fluxes at a nominal wavelength (in this case, A5100A) as the main convergence criterion. Eor the calculation 
of an individual spectrum, the radial Thomson depth the radius of the Comptonizing region r^, the disk inner radius r^ain, the central 
black hole mass, the accretion rate, and initial estimates of the electron temperature and X-ray luminosity are taken from our semi- 
analytic calculations of §3. Eor the first iteration, we assume a nominal value for the disk albedo, and we use the expressions from 
Zdziarski et al. (1999), modified to include the contribution from local viscous dissipation, to determine the effective temperature as a 
function of disk radius. On subsequent iterations, the model and observed optical fluxes are compared, and the electron temperature 
(and therefore the X-ray luminosity) is adjusted. Since we are computing the Compton reflection spectrum in the disk at each 
iteration, we obtain the disk albedo as a function of radius which we use as input for the subsequent iteration. Convergence has been 
achieved when both the model optical flux agrees with the observed value and the X-ray flux has converged to within a specified 
tolerance—typically 5% for each. As with our semi-analytic calculations, our Monte Carlo calculations do not consider pair balance 
or the presence of a non-thermal component in determining the properties of the Comptonizing plasma. 
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M dM/dt 

Fig. A1.— Equilibrium electron temperature, ksTg, versus MM for the 1989-90 epoch 2 data, plotted as different clusters of points from the upper left to the 
bottom right. Each cluster corresponds to a different chosen value of MM = (0.5,1.5,2.5,3.5,4.5,5.5,6.5,7.5) x 10^ yr“^. Within each cluster the different 

point symbols correspond to different black hole masses of M-j = A//( 1 O^M 0 ) = (1,2,3,4,5,6,7), as indicated in the figure. Note that ksTe depends primarily 
on MM, with only a weak dependence on M arising because of the resulting variation in rj (see eq. [2]). The solid curve shows the relation between kgTe and 
MM after we have averaged over this weak mass dependence within each cluster of points. The dashed curves represent constraints from global energy balance, 
for different values of black hole mass as indicated by the point symbols at the ends of these curves. Each of these dashed curves shows kgTe vs MM^nin, where 
Mmin = iLx+Lir)/ir)c^) and 77 = 0.5. For a given black hole mass, the intersection of the dashed curve with the solid curve gives the minimum accretion rate and 
maximum electron temperature (see §4) that are consistent with equations 15a-e. 





12 



M dM/dt 


Fig. A2.— Electron temperature versus MM for all nine 1989-90 epochs. The kgTe values for a given MM have been averaged over the seven black hole masses 
considered in Figure Al. The kgTe-MM curves for each epoch are identified by the combination of linestyles and symbols as shown in right-hand side of this figure. 
The dotted vertical lines indicate the minimum value of MM for each mass for the epoch 1 data, which are the most constraining. These lines correspond to black 
hole masses M = (1,2,3,4.5,6.7) x 10^ Mq, ordered left-to-right and labeled near the top of the plot. The dotted horizontal lines show our assumed limits on the 
electron equilibrium temperatures, - 30 keV and ksTg^.^^ = mgC^ = 511 keV. If we require ksTg > 30 keV for all nine epochs, then the ksTg-MM curves for 

epochs 7 and 9 constrain the black hole mass to be M < 2 x 10^ Mq and MM < 2.5 x 10^ MQ^yr“^. The arrows near the top of the plot indicate the minimum 
Mlf-values averaged over the constraints for all nine epochs and correspond to M = (1,2,3,4,5,6,7) x 10^ Mq, ordered left-to-right. Therefore, taken as a whole, 
the opticayiUE/Ginga data require M < 3 x 10^ Mq. 
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Fig. A3.— Spectral energy distributions which have been “fit” to the optical, UV, and X-ray data for the nine 1989-90 epochs. The optical/UV data are plotted 
as the open diamonds and the power-law fits to the 2-10 keV Ginga data are plotted as the thick line segments. The solid curves are the SEDs obtained from our 
semi-analytic calculations described in §3, and the histograms are the thermal Comptonization/thermal reprocessing Monte Carlo results using the parameters, r.y, 
r^in. and Tr = T, from the semi-analytic calculations (see §4). All these SEDs have been computed assuming M = 1 x 10^ Mq and M = 0.02 Mq yr“^ 












14 



Fig. A4.— ksTe versus rmin for the seven epochs of 1998 ground-based/RXTE data. The mass and accretion rate have been fixed at M = 1 x 10^ Mq and 
M = 0.02 Mq yr“^. The horizontal dotted lines are our assumed upper and lower limits for the electron temperature, ksTef^^tx = 511 keV and - 30 keV. The 

individual curves have been obtained using the 1999 RXTEf^CA response matrices. The lower limit the epoch 2.2 curve imply a maximum disk inner 

radius of < 2.4 x 10^“^ cm (leftmost vertical dotted line). By contrast, the upper limit the epoch 3.2 curve require a minimum disk inner radius, 

^min ^ 3.2 X 10^^ cm (rightmost vertical dotted line). Taken together, these constraints imply that the disk inner radius had to vary significantly over the course of 
the 1998 observations. The hashed region shows the range of ksTe values at each r^jn that are obtained using the 2002 RXTE/PCA response matrices. For these 
fits, a single value of in the range (3.i-4.7)x 10^“^ cm may be taken for all seven epochs and will be consistent with our nominal upper and lower electron 
temperature limits. 
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Fig. A5.— Model spectral energy distributions for the 1998 ground-based/RXTE data. Using the 1999 response matrices, SEDs have been computed for two 
inner disk radii, rmin = 2.4 x 10^^ cm (solid curves) and rmin = 3.2 x 10^^ cm (dashed curves). The differences in thermal Comptonization cut-off energy and the 
location and magnitude of the spectral peak in the UV band can be clearly seen. In epochs 2.2-2.4, the equilibrium electron temperatures are all < 30 keV if one 
assumes a disk inner radius of = 3.2 x 10^“^ cm. For an exponentially cut-off power-law with Ec = IhgTe, these spectra differ substantially from the single 
power-law component inferred for the 2-10 keV band from the RXTEf^CA data (thick solid line segment). The SED calculations using the 2002 response matrices 
are shown as dot-dashed curves and a thinner line segment is plotted showing the 2-10 keV power-law for these data. 













16 



A-1 


Fig. A6 .— Photon spectral index versus Compton amplification factor minus unity. The filled triangles are from the Monte Carlo calculations shown in Figure A3 
that were computed using the parameters for the 1989-90 data. The other symbols correspond to separate sets of Monte Carlo calculations for geometries with 
values of rs/r^m. = 0.5, 0.75, 1, 1.25, and 1.5, as indicated in the figure. Each set of Monte Carlo calculations had M = 2 x 10^ Mq and different values of M = 0, 
0.01, 0.02, 0.04, and 0.08 Mq yr“^, with larger M values corresponding to larger values of spectral index F. (For larger values of M, smaller fractions of soft 
photons from thermal reprocessing are required, and so softer X-ray spectra are produced.) The solid curve is the parameterization found by Malzac et al. (2001) 
which has Fq = 2.15 and <5 = 14 (see eq. 18]), while the dashed curve is the best fit to the Monte Carlo data yielding Fq = 2.06 and 5 = 15.6. 
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Fig. A7.— The ratio of the luminosity of thermally reprocessed seed photons, L 5 ,repr, to the apparent X-ray luminosity, L;c,app, versus the ratio of Comptonizing 
sphere radius, rs, to disk inner radius, rmin- The data points are the Monte Carlo calculations, and the various symbols correspond to the same datasets, rs/r^in 
and Ginga data, as in Figure A6. The solid curve is the relation obtained assuming the thermal Compton emission is isotropic. The dashed curve is the solid curve 
multiplied by ^ = 1.5. 
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M dM/dt (MgV^^) 

Fig. A8 .—ksTe versus MM for two different observer inclinations, cos/ = 0.866 (upper) and cos / = 0.5 (lower). These plots illustrate the sensitivity of the limits on 
the black hole mass and accretion rates to the assumed disk inclination. For cos / = 0.866, the dotted vertical lines are the minimum MM-values for black hole masses 
M = 0.2,0.4,0.6,0.8,1.2,3 x 10^ Mq, ordered left-to-right, while for cos / = 0.5, the dotted vertical lines correspond to MM limits for M = 1,2,3,4,5 x 10^ Mq. 
In both cases, these limiting values are from the epoch 1 curve. The horizontal dotted lines indicate = 30 keV and kgTe^.^ = 511 keV. 
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Table A1 

Optical and UV continuum fluxes and X-ray continuum parameters from the 1989-90 GROVND-BASED/IUE/GingaDATA. 


epoch 

JD-2440000 

fsioo" 

±2670^ 

±1840" 

±1350" 

± 210 '’ 


pc 

RC,d 

1 

7535-7539 

0.78 ± 0.08 

2.22 ± 0.19 

4.52 ± 0.22 

7.45 ± 0.33 

3.39 

1.56 

±0.19 


2 

7556 

0.62 ± 0.06 

2.09 ± 0.20 

3.91 ±0.22 

6.18 ±0.33 

5.43 

1.72 

±0.03 

0.58 ±0.20 

3 

7685 

0.67 ± 0.06 

1.72 ±0.23 

2.88 ±0.31 

3.94 ± 0.46 

4.71 

1.71 

±0.10 

0.44!||:« 

0.28i24 

4 

8036-8037 

0.50 ± 0.06 

1.31 ±0.36 

2.56 ± 0.30 

3.84 ± 0.35 

3.87 

1.68 

±0.06 

5 

8039-8041 

0.51 ±0.07 

1.44 ±0.25 

2.85 ± 0.36 

4.25 ± 0.42 

3.53 

1.78 

±0.06 

o-9mi 

6 

8047 

0.49 ± 0.07 

1.40 ±0.19 

2.59 ± 0.30 

4.21 ± 0.35 

2.94 

1.75 

±0.08 

0.97!«:g 

7 

8056 

0.50 ± 0.06 

1.53 ±0.17 

2.80 ± 0.21 

4.35 ± 0.24 

4.95 

1.74 

±0.04 


8 

8068 

0.39 ± 0.05 

1.10±0.13 

1.64 ±0.30 

2.23 ± 0.35 

2.44 

1.53 

±0.04 

<0.27 

9 

8077 

0.34 ± 0.05 

0.96 ± 0.20 

1.30 ±0.30 

1.74 ±0.35 

1.93 

1.71 

±0.14 

0.87!“;^^, 


"Units are 10 erg cm ^A The flux uncertainties are l-cr. 

^Units are 10“^^ erg cm“^s“^. Throughout this paper, the 2-10 keV flux is computed for the underlying X-ray continuum which is assumed to be a cut-off power-law 
in the PEXRAV model. Therefore, these fluxes exclude the Fe Ko; emission, the Compton reflection component, and corrects for the effects of absorption. 
‘^Uncertainties are computed for Ax^ = 2.706. 

‘^Compton reflection fraction. 


TableA2 

Model parameters lit to the 1989-90 GROUND-BASED//([/£/GmgflDATA ASSUMING M= 1 X lo’ Mq, M = 0.02 Mq yr“', and 

cos; = 0.7. Note that rmm and To depend only on cos;. 


epoch 

Tmin" 

To (10" K) 

r ^ 

ksTe (keV) 

r 


T 

^tr 

1 

2.20 

4.0 

2.07 

683 

0.13 

4.65 

0.11 

2 

1.85 

4.2 

2.07 

153 

0.72 

3.03 

0.13 

3 

3.36 

2.8 

3.70 

66 

1.64 

2.10 

0.07 

4 

2.31 

3.3 

2.42 

76 

1.54 

1.84 

0.10 

5 

2.07 

3.6 

2.47 

96 

1.05 

1.68 

0.11 

6 

1.97 

3.7 

2.22 

226 

0.44 

1.77 

0.12 

7 

1.96 

3.7 

2.19 

37 

2.46 

1.90 

0.12 

8 

2.57 

2.8 

2.17 

85 

1.82 

1.44 

0.09 

9 

2.62 

2.7 

2.76 

50 

2.05 

0.80 

0.09 


“Units are 10*'* cm. 
*’Units are 10'*'* erg s“*. 
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Table A3 

Optical continuum fluxes and X-ray continuum parameters from the 1998 ground-based/RXTE data. 


epoch JD-2450000 


Fsioo^ 1999 PCA responses” 

77 d “pe DC 

^210 ^ ^ 


2002 PCA responses'^ 

T7 d pe oe 

^210 i ^ 


1.0 

979.96- 

-980.40 

1.16 

+ 0.06 

8.0 

1.86 + 0.03 

0.34!!J;;« 

7.9 

1.72 + 0.03 

0.20 

2.1 

985.09- 

-985.40 

1.17 

+ 0.09 

7.8 

1.86 + 0.03 

0.38tj^ 

7.7 

1.72 + 0.03 

0.24 

2.2 

985.69- 

-986.41 

1.18 

+ 0.09 

9.8 

1.93 + 0.02 

0.46 + 0.12 

9.5 

1.80 + 0.02 

0.33 

2.3 

986.65- 

-987.40 

1.16 

+ 0.09 

9.5 

1.90!«:02 

0 37+0.10 
'-0.04 

9.3 


0.25 

2.4 

3.1 

987.63- 

994.36- 

-987.86 

-994.89 

1.24 + 0.09 
1.03 + 0.04 

9.1 

6.9 

1 oa+ 0.05 
l-oU-0 04 

0 46+0-18 

0.29ii 

8.8 

6.8 

1 ^q+0.05 
^•^^-0.04 

I •DD_o 04 

0.34 

0.17 

3.2 

995.64- 

-996.13 

1.01 

+ 0.08 

5.6 

1 oa+0.03 
^•^^-0.02 

0.38ti« 

5.5 

1 ^c+0.04 
^•^^-0.03 

0.23 


TOTT 

11? 

-0.14 

+0.11 

- 0.10 

+0.11 

- 0.10 

+0.28 

- 0.20 

+0.24 

-0.17 

+0.17 

0.14 


"Units are 10“^"^ erg l-cr uncertainties. 

^Using PCARSP V2.36 and PCARMF v3.5. 

•^Using PCARSP v8.0 and PCARMF v8.0. 

^Units are 10“^^ erg cm^s“^. 

^Uncertainties are computed for Ax^ = 2.706. 


Table A4 

Model parameters fit to the 1998 ground-based/RXTE data for rmin = 2.4 x 10^"^ cm and r^m = 3.2 x 10^"^ cm (in 

PARENTHESES). THESE PARAMETERS WERE DERIVED FROM SPECTRAL FITS THAT USED PCA RESPONSE MATRICES GENERATED BY 

PCARSP V2.36 AND PCARMF v3.5 (Chiang ET AL. 2000). 


epoch 

Vs (10‘'* cm) 

keTe (keV) 

T 

L/ 

1.0 

3.78 (5.02) 

554(136) 

0.09 (0.63) 

5.15 (4.04) 

2.1 

3.79 (5.02) 

684(168) 

0.06 (0.49) 

5.21 (4.09) 

2.2 

4.64 (5.99) 

27 ( 12) 

2.26 (3.80) 

4.28 (3.94) 

2.3 

4.24 (5.52) 

89 ( 12) 

0.90 (3.99) 

4.57 (3.75) 

2.4 

4.38(5.81) 

204 (27) 

0.36 (2.31) 

4.92 (3.86) 

3.1 

3.24(4.31) 

676 (229) 

0.07 (0.37) 

5.00 (3.95) 

3.2 

3.24(4.31) 

1552 (526) 

0.02 (0.11) 

4.87 (3.84) 


"Units are 10^ erg s ^ 


Table A5 

Model parameters fit to the 1998 ground-based/RXTE data for rmin = 3.2 x lO'"^ cm. These parameters were derived 
FROM SPECTRAL FITS THAT USED PCA RESPONSE MATRICES GENERATED BY PCARSP V8.0 AND PCARMF V8.0. 


epoch 

o 

o 

keTe (keV) 

T 

L/ 

1.0 

3.72 

346 

0.25 

5.57 

2.1 

3.72 

393 

0.21 

5.63 

2.2 

4.34 

135 

0.69 

4.85 

2.3 

4.07 

164 

0.59 

5.05 

2.4 

4.16 

305 

0.26 

5.48 

3.1 

3.43 

285 

0.38 

5.05 

3.2 

3.39 

471 

0.19 

4.98 


'Units are 10^ erg s ^ 



